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I .  INTRODUCTION 

During  the  past  several  years  yawsondes  have  been  used  to  measure 
in-flight  dynamical  behavior  of  projectiles  over  long  flight  paths1  . 

The  yawsonde  produces  large  volumes  of  data  on  the  solar  aspect  angle 
and  spin  of  the  shell  during  flight.  Yawsonde  data  are  two-dimensional. 
The  angle  between  the  projectile  axis  and  the  sun  as  well  as  the  spin 
are  determined  along  the  flight  path.  The  actual  angular  motion  of  the 
shell,  however,  is  three-dimensional  consisting  of  pitch,  yaw,  and  spin. 
The  basic  problem  is  to  extract  the  three-dimensional  behavior  of  the 
shell  from  two-dimensional  data.  Additional  information  on  the  tra¬ 
jectory  and  meteorological  conditions  are  required  as  inputs. 

The  conventional  method  for  reducing  data  from  free  flight  tests 
is  to  fit  solutions  to  the  equations  cf  motion  to  the  range  data2. 

Such  procedures  assume  that  the  equations  of  motion  can  be  solved  in 
closed  form  or  in  an  approximate  way  (quasi -1 inear) .  The  conventional 
method  can  be  applied  to  yawsonde  data  using  a  technique  called  WOBBLE3. 

A  modified  version  of  WOBBLE  has  been  used  with  limited  success1.  Still 
another  method,  developed  by  Haseltine4,  uses  an  analog  computer  set  up 
with  the  equations  of  motion.  The  operator  adjusts  the  coefficients  of 
the  equations  until  the  motion  predicted  by  the  analog  computer  agrees 
with  the  flight  data. 

Chapman  and  Kirk5  recently  developed  a  method  for  fitting  differen¬ 
tial  equations  of  motion  directly  to  the  data  using  a  unique  approach 
for  adjusting  the  coefficients  of  the  equations.  This  method  does  not 
require  that  the  equations  of  motion  have  an  analytical  solution.  Whyte 
and  Beliveau  have  investigated  the  method  of  Chapman  and  Kirk  and  applied 
it  to  data  from  aerodynamic  spark  range  firings  with  excel lecnt  results6. 

In  this  report  it  will  be  shown  that  the  numerical  integration 
technique  can  be  applied  to  yawsonde  data  to  obtain  linear  aerodynamic 
coefficients  over  extended  portions  of  the  flight  of  a  single  projectile. 
Actual  reductions  were  performed  on  a  variety  of  projectile  types.  Here 
the  discussion  is  limited  to  the  work  done  on  the  155mm,  M549  rocket 
assisted  projectile.  This  projectile  was  chosen  so  that  direct  compar¬ 
isons  could  be  made  with  the  results  of  aerodynamic  range  and  wind 
tunnel  tests. 

This  report  will  discuss  the  aerodynamic  equations  for  angular  and 
swerving  motion  and  the  numerical  integration  method  itself.  The  re¬ 
sults  of  axial  force  and  angular  motion  reduction  will  be  presented  and 
compared  with  linear  theory  reduction  of  range  an  wind  tunnel  tests. 

The  effete  of  including  wind  and  velocity  dependent  coefficients  in  the 
method  is  described  and  discussed. 


^References  are  listed  on  page  29. 
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II.  THEORETICAL  CONSIDERATIONS 


Linearized  theory  has  been  used  to  describe  the  mo"  ion  of  symmetric 
missil  with  outstanding  success  during  the  past  two  decades.  The 
theory  xs,  however,  restricted  to  small  angle  motion  with  small  velocity 
and  spin  changes  and  nearly  linear  aerodynamic  coefficients.  Nonlinear 
coefficients  may  be  determined  by  linear  theory  from  overlapping  fits  if 
sufficient  data  are  available  or  if  sufficient  number  of  range  tests  of 
the  same  projectile  are  made. 

None  of  the  above  restrictions  apply  to  the  method  of  Chapman  and 
Kirk.  This  method  is  conceptually  and  practically  attractive  since  aero¬ 
dynamic  coefficients  can  be  determined  by  fitting  the  differential  equa¬ 
tions  of  motion  directly  to  the  observed  angular  motion  data.  The  force 
coefficients  can  be  obtained  by  fitting  the  equations  of  swerving  motion 
directly  to  positional  data.  The  general  method  can  be  summarized  briefly 
as  follows.  The  equations  of  motion  appropriate  to  the  behavior  of  the 
projectile  are  numerically  integrated  using  estimated  values  of  the  co¬ 
efficients  and  the  initial  conditions.  The  same  equations  of  motion  are 
partially  differentiated  on  each  coefficient  forming  a  set  of  parametric 
equations.  The  parametric  equations  are  then  numerically  integrated  to 
obtain  values  of  the  partial  derivatives  with  respect  to  each  coefficient. 
A  differential  corrections  equation  is  set  up  from  a  Taylor  expansion  of 
the  dependent  variables  of  the  equations  of  motion.  The  data  are  then 
compared  to  the  computed  values  of  the  dependent  variables  in  a  least 
squares  sense  and  corrections  to  the  coefficients  are  obtained  from  the 
differential  corrections  equations.  The  coefficients  of  the  equations 
of  motion  are  adjusted  by  these  corrections  and  the  process  begins  again 
with  new  values  of  coefficients.  The  method  iterates  until  convergence 
is  achieved. 

A.  Equations  of  Motion 

The  equations  of  motion  of  a  symmetric  missile  have  been  derived  in 
a  fixed  plane  coordinate  system7.  Hie  rotation  of  the  earth  is  ignored 
and  the  aerodynamic  coefficients  are  assumed  to  have  the  form  of  a  power 
series  (as  functions  ot  sin  a). 

a  +  a  e2  +  a  e4  +  a  (V  -  V  )  (1) 

0  3  5  V  o 


b  +  b  e2  +  b  c4  (2) 

0  3  5 


d  +  d  e2  (5) 

0  2 
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The  static  and  Magnus  coefficient  expansions  include  fifth  order  co¬ 
efficients  and  may  even,  in  some  cases,  require  higher  order  terms. 

The  motion  of  the  center  of  mass  of  the  projectile  is  described  by  the 
force  equations  (4-6) . 

u  =  -  q  A  Cx/m  +  g  sin  %  -  ®fP  w  ♦  % cos  % v  (4) 


l£.C«l— — 
m  a  V  m  2V  xp  V 
w  w  r  a  w 


%  cos  e£p  [“  +  w  tan  V 


w  =  -q  AC.  Ww  -  q  A  pL  C,  Vw  -  g  cos  6 . 

m  a  V  m  2V  p  V 
w  w  *a  w 

.  ifp  v  sin  0fp  *  8^  u  (6) 

The  angular  motion  of  the  projectile  is  described  by  the  moment  equations 
(7-9). 

%  *  *  iep  sin  efp  .  ifp  8fp  cos  9fp  (7) 

x  ^Vw  * 


where  the  spin  is  p  =  ^  sin  0^ 


%  ■  « 


AL  fa  +  a  e2  +  a  e4  +  a  (V  -  V  )1  w  JV 

-  1“  3  5  v  "J 

+  qAL2  f  d  +  d  e2 1  0  - 

zv  L 0  2  J  * 

♦  q  AL^  ft  .be2,  b  c"l  \ 

2V  I  I  0  3  5  J  V 

w  «•  J  w 

-  P  l*  cos  0*  -  *2  cos  0-  sin  0- 


d 


i 


I 


wt 


"fp 


[-  <t  “  |«0  *  %  ^  -  %  ‘4  ♦  %  <v  -  y]  vv. 


-fp  cos  % 


♦  q  AL2  fd  +  d  e2  1  4, 

iy  l*  2  J 

+  q  AL^  Tb  +  b  e2  +  b  Ww 

2VI  |  3  3  5  J  V 

w  -  J  w 

*  2  5fp  *fp  sin  8fp  *  %  t  iij/cos  6^ 


(9) 


0fp  and  -4»fp  are  angles  which  define  missile  orientation  in  the  earth 


fixed  plane  coordinate  system.  0^  describes  a  rotation  about  the 


missile  y-axis  and  ^  is  a  rotation  about  the  earth  z-axis.  It  is 


convenient  to  express  the  force  equations  (4-6)  in  an  earth  fixe;!  co¬ 
ordinate  system.  The  transformations  are  equations  (10-12) . 


Vx  =  u  cos  0fp  cos  ^  -  v  sin  +  w  sin  0fp  cos  if»fp 


Vye  =  u  cos  efp  sin  ♦fp  +  v  cos  ^  +  w  sin  0^  sin  ^ 


(10) 

(ID 


V  =  -  u  sin 
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fp 
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this  point  we  wish 

to  take 

the  effect  of  wind  into  account. 
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Differentiation  of  the  transformation  equations  (10-15)  and  substitution 
into  the  force  equations  (4-6)  gives  the  equations  of  motion  of  the 
center  of  mass  in  an  earth  fixed  system: 


■  -  q 


\  c2l cos  % cos  % 


A  fC  +  C  e2l 

5  [  XD  x2  J 
*  ’i1*.  (£sin  '  r sin  6fp cos  *fp) 

-  qA  pL  C.  \  sin  ip.  +  vw  sin  0„  cos 

-  2\  \  [\  fp  \  *  rpj 

h  [\  (v“  -  v] 


COS  9_  COS  llir 

fp  Yfp 


(16) 


Cx2  e2  |  cos  6fp  sin  *fp 


=  q 


5S  [VV‘] 

*  <1  AC.  f Vw  cos  ,(»-  +  ww  sin  0-.  sin  ip.  \ 

»  \JT  fp  \T  fP  fP| 

+  q  A  pL  C,  ( Ww  cos  *_ 

r  ^ 

w  *a  \  w 

-  qA  fc  (V  -  V  )1 

»  [  \  w  0  J 

5  ['*.  *  S 1 


-  w  sin  6  Sin  *_ 

V~  fp  *P 

w 


cos  e.  sin  ip. 
fp  fp 


■) 


(17) 


sm  0 


q  a  _w  cos  e . 

-  fp 


q  A  pL  C,  w  cos  0  _ 

IT  2V“  S  V"  fp 
v*  a  w 


(18) 
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where 


V  =  /v2  ♦  V2  +  V2  \  ^ 

\  Xe(w)  VeM  *e(w)/ 


,  and  q  =  *s  p  Vj 


axial  force  coefficient  has  been  expanded  as 
C  -  C  +  C  e2 

XXX 
0  2 


and  the  normal  and  Magnus  force  coefficients  can  be  expanded 

C.,  =  n  +  n  e2 

X  0  3 


C,  =  m  ♦  m 

V  0  3 

a 


A  Mach  number  dependent  axial  force  coefficient,  ,  has  been  included 

v 

in  equations  (16-18)  since  large  velocity  changes  are  cften  encountered 
in  flight  data. 

Time  is  the  independent  variable  in  the  above  equations  rather 
than  distance  along  the  traj  cctory  (arc  length) .  Distance  is  a  more 
natural  base  for  data  taken  in  aerodynamic  spark  ranges,  while  time  is 
the  natural  variable  for  yawsonde  and  radar  data.  The  method  of  Chapman 
and  kirk  is  now  used  in  a  computer  program  called  IIEEVE  to  best  fit 
equations  (16-18)  to  the  radar  position-time  data.  Estimated  values  of 
6  and  v  are  used  to  obtain  fitted  values  of  position,  velocity,  axial 
force,  normal  force,  and  Magnus  force.  Before  the  angular  motion  equa¬ 
tions  (7-9)  can  be  fitted  to  the  yawsonde  data,  a  relationship  between 
the  solar  aspect  angle  and  the  angles  6^,  must  be  derived. 

B.  Solar  Aspect  Angle  Relation 

Yawsonde  data  give  one-dimensional  angular  motion  of  the  missile 
with  respect  to  the  sun.  The  solar  aspect  angle,  a,  is  the  angle  be¬ 
tween  the  longitudinal  axis  of  the  projectile  and  a  vector  directed 

from  the  projectile  to  the  sun.  Let  t  ,  t.  ,  t  be  the  direction  cosines 

12  3 

of  the  solar  vector  in  an  earth  fixed  coordinate  system  (Fig.  1).  e^, 

and  locate  the  axis  of  the  projectile  in  the  same  coordinate  system. 

The  solar  aspect  angle  is  obtained  from  the  dot  product  of  a  unit  solar 
vector  and  a  unit  vector  along  the  missile  axis. 


i 


t 

i 


i 


j 


i 


cos  a  =  t  cos  6.  cos  1 l>-  +  Z  cos  9,  sin  +  t  sin  6_  (22) 

1  fp  fp  2  fP  3  fP 

Equation  (22)  provides  the  required  relationship  between  the  sola?  aspect 
angle,  the  orientation  of  the  sun,  and  the  angles  of  the  missile  axis 
relative  to  the  earth. 

The  angular  motion  of  the  projectile  is  fitted  to  the  angular  equa¬ 
tions  (7-9)  in  a  least  squares  sense  by  comparing  computed  values  of  a 
using  equation  (22)  to  the  measured  values  of  a.  The  method  of  Chapman 
and  Kirk  is  used  to  adjust  the  coefficients  of  the  angular  motion  equa¬ 
tions  until  the  probable  error  of  fit  is  a  minimum. 


III.  YAWSONDE  REDUCTION  PROCEDURE 

In  order  to  obtain  aerodynamic  coefficients  from  yawsonde  flight 
data,  the  center  of  mass  met ion  of  the  projectile  must  first  be  well 
determined.  Measurements  of  the  physical  properties  of  the  shell  are 
taken  before  firing  and  meteorological  data  on  density,  temperature , 
and  wind  velocities  are  noted  at  the  time  of  firing.  First  a  reduction 
on  the  position  time  data  is  made  to  obtain  a  best  description  of  the 
trajectory.  The  results  of  this  reduction  are  used  as  inputs  to  the 
angular  motion  reduction.  If  necessary,  the  results  of  the  angular 
motion  reduction  could  be  fed  back  into  the  position-time  reduction 
for  improved  accuracy. 

The  sequence  of  data  reduction  using  the  numerical  integration 
technique  Is  as  follows.  First  the  raw  radar  position-time  data  are 
reduced  to  obtain  a  drag  profile.  The  inputs  to  this  swerve  reduction 
are  measured  values  of  position  (xe,ye,z£)  and  time,  physical  properties 

of  the  projectile,  meteorological  data  on  winds,  density  and  temperature, 
and  the  direction  cosines  of  the  solar  vector  in  an  earth  fixed  co¬ 
ordinate  system.  The  swerve  reduction  also  requires  estimated  values 
of  axial  force,  normal  force,  Magnus  force,  spin  profile,  spin  damping, 
and  static  moment  coefficient.  Such  estimates  can  be  obtained  from  the 
known  aerodynamics  of  similar  projectiles  or  from  ballistic  range  measure¬ 
ments.  The  swerve  reduction  provides  an  output  of  axial  force  (related 
to  drag)  versus  Mach  number,  the  spin  damping,  and  best  fit  computed 
values  of  position,  velocity,  dynamic  pressure,  and  spin,  all  as  func¬ 
tions  of  time. 

The  results  of  the  swerve  reduction  are  used  as  inputs  to  the 
angular  reduction  along  with  the  solar  aspect  angle  data  from  the  yaw¬ 
sonde.  Hie  angular  reduction  also  requires  the  physical  properties  of 
the  projectile  and  the  meteorological  data  as  inputs.  Initial  esti¬ 
mated  values  of  the  static  moment,  Magnus  moment,  and  damping  moments 
must  also  be  furnished  The  angular  reduction  provides  an  output  of 
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best,  fitted  values  of  solar  aspect  angle  versus  time  plus  the  aero¬ 
dynamic  coefficients  for  static,  Magnus  and  damping  moments.  Since  the 
yawsonde  data  usually  span  extensive  portions  of  the  trajectory  includ¬ 
ing  many  periods  of  precession  and  nutation,  the  data  arc  fit  section- 
ally.  The  results  of  sectional  fitting  are  average  valuer,  of  the  aero¬ 
dynamic  coefficients  over  the  data  intervals.  The  entire  reduction 
procedure  is  summarized  in  the  following  table. 


Table  I.  Reduction  Inputs  and  Outputs 


Swerve  Reduction 

Angular  Reduction 

Inputs 

Measured  x  ,y  ,z  ,  t,  mass, 
e  ye  e 

diameter,  air  density,  wind 
velocity,  direction  cosines 
of  sun,  spin  profile. 

Estimated  C^C^Cy  ,C£  ;Cffl 

P  P 

Direction  cosines  t  ,t  ,t  . 

12  3 

Measured  solar  aspect  angle, 
mass,  moments  of  inertia, 
air  density,  wind  velocity. 
Fitted  p,  q,  V,  Vw- 

Estimated  C^Cy  ,Cm,Cn  ,Cffl 

P  p  q 

Outputs 

Computed  xe,ye,ze,V,Vw,q,p. 
Fitted  Cx,  C£ 

Computed  a  versus  t. 

Fitted  ,C  . 

P 

q  P 

IV.  RESULTS  AND  DISCUSSION 

The  results  of  fitting  for  the  translational  avid  angular  motion  of 
a  155mm,  M549  rocket  assisted  projectile  ore  discussed  below.  This 
projectile  was  fired  at  Wallops  Island,  Virginia,  cr.  March  11,  1971,  at 
2009  hours  Universal  Time  from  a  1/20  twist  tube.  The  muzzle  velocity 
as  determined  by  smear  cameras  was  321  aeters/sec.  The  gun  elevation 
was  490  mils.  The  physical  characteristics  of  the  shell  and  the  wind 
data  are  summarized  as  follows: 

Weight  =  42.75  kg 

Location  of  c.g.  =  32.21  cm  from  base 

Axial  moment  of 

inertia  -  0.1477  kg-m2 

Transverse  moment  of 

inertia  =  1.8009  kg-m2 


Diameter 


15.5  cm 


i 


Length 

Mean  wind  speed 


=  87.4  cm 

=  12.8  ra/sec 


Mean  wind  direction  =  300  d_*g.  from  true  North 

An  outline  oc  the  projectile  is  shown  in  Figure  2. 

A.  Swerve  Reduction 

A  swerve  reduction  was  performed  on  the  data  obtained  from  an  FPS-J.6 
radar  which  tracked  the  projectile  from  5.3  seconds  after  launch  to  im¬ 
pact.  In  the  swerve  reduction  the  yaw  of  repose  was  estimated.  Swerve 
reductions  were  performed  both  with  and  without  wind  effects.  The  re¬ 
sults  of  the  swerve  reductions  are  given  in  Table  II.  The  probable 
error  in  each  parameter  is  shown  in  parentheses  below  the  parameter. 


4 4 


Table  II.  Results  of  Swerve  Reduction 


Muzzle 

m/s 


Impact 

ra/s 


Probable 

Error 


t: 


constant 


The  muzzle  velocity  without  wind  effects  agrees  well  with  the  smear 
camera  value  of  muzzle  velocity.  The  errors  in  the  fit  decrease  when 
wind  effects  are  taken  into  account.  In  order  to  accomplish  the  swerve 
reduction,  values  of  normal  force,  Magnus  force,  and  nonlinear  axial 
force  coefficients  were  estimated  and  held  co<istant  during  the  reduction. 

The  spin  damping  coefficient-  ,  was  determined  for  tne  flight 

P 

by  matching  a  computed  spin  history  to  a  flight  spin  history'  derived 
from  the  yawsonde  data.  In  this  manner,  C£  was  found  to  be  -  0.025 

P 

±  .001. 


■A 


f  1, 


B.  Angular  Reduction 

A  portion  of  the  data  from  the  yawsonde  on  board  the  projectile  is 
shown  in  Figure  3-  The  data  are  solar  aspect  angle  versus  time  between 
5.0  and  17.0  seconds  of  the  flight.  The  data  clearly  show  a  fast  fre¬ 
quency  and  a  slow  frequency  which  can  be  identified  as  nutation  and 
precession.  The  nutational  amplitude  damps  out  rather  quickly  while 
the  precessional  amplitude  grows  slowly.  In  order  to  obtain  a  reason¬ 
able  value  of  damping  moment  at  least  several  cycles  of  the  precession 
must  be  included  in  the  fit.  The  static  moment  will  be  the  dominant 
factor. 

A  series  of  reductions  were  performed  with  variations  in  section 
length  and  data  density.  Data  density  refers  to  the  number  of  data 
points  per  unit  time  included  in  the  fit.  Reductions  were  done  with 
and  without  wind  effects.  Various  combinations  of  linear  ano  nonlinear 
coefficients  were  tried.  The  results  of  some  reduction  attempts  are 
listed  in  Table  III.  On  some  reductions  a  velocity  dependent  static 
moment  coefficient,  ay,  was  sought.  On  other  reductions,  the  ratio  of 

moments  of  inertia  was  a  parameter  of  the  fit. 

The  coefficients  were  not  substantially  affected  by  changes  in 
section  length  or  data  density  as  can  be  seen  by  comparing  runs  2,  3. 
and  7.  The  maximum  length  used  was  14.41  seconds  (run  7).  For  this 
length  the  probable  error  of  fit  was  0.31  degree.  Runs  8-10  repre¬ 
sent  an  attempt  to  extract  nonlinear  coefficients.  A  comparison  be¬ 
tween  runs  1  and  7  shows  the  effect  of  wind.  A  change  of  10%  in  static 
moment  is  observed  and  is  considered  significant.  The  changes  in  damp¬ 
ing  and  Magnus  moments  are  small  and  not  significant. 

The  average  total  angle  of  attack  for  this  flight  was  computed  to 
be  i.72  degrees  which  represents  small  amplitude  angular  motion.  Runs 
2-7  show  the  linear  static  moment  coefficient  varying  from  4.64  to 
4.76  at  an  average  Mach  number  of  .74  over  this  part  cf  the  flight. 

The  fit  for  run  7  can  be  seen  in  Figure  4  where  the  computed  and  mea¬ 
sured  solar  aspect  angle  are  plotted  against  time  of  flight.  For  the 
sake  of  short  figures,  only  a  section  of  the  fit  interval  is  displayed. 

In  Figure  5  the  values  of  static  moment  coefficients  from  run  7 
(linear)  and  run  10  (nonlinear)  are  compared  to  data  from  wind  tunnel 
tests  performed  on  the  M549  projectile.  The  total  static  moment  at  an 
average  angle  of  attack  of  1.72  degrees  is  computed  for  an  average 
Mach  number  of  0.74.  Both  the  linear  fit  moment  and  the  nonlinear  fit 
moment  agree  well  with  wind  tunnel  results.  To  show  the  effect  of  not 
including  wind,  the  static  moment  result  of  run  1  is  also  shown.  The 
point  is  noticeably  away  from  the  wind  tunnel  data. 

Wind  tunnel  Magnus  moment  measurements  at  the  above  values  of  Mach 
number  and  angle  of  attack  give  a  C  =  .0125  while  the  Magnus  moment 
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Derived  from  MS49  Yawsonde  Data 
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Table  III.  Aerodynamic  Coefficients  Derived  from  M549  Yawsonde  Data  (continued) 
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computed  from  run  7  is  -  .0124.  The  damping  moment  coefficient  from 
wind  tunnel  tests  is  -  5.1  while  run  7  shows  +3.4.  In  light  of  the 
yawsonde  data  shown  in  Figure  3,  it  would  seem  that  Magnus  and  damping 
moments  cannot  be  readily  extracted  from  small  amplitude  motion  where 
the  nutational  mode  dominates  the  data. 


CONCLUSIONS 


A  technique  of  numerical  integration  has  been  developed  and  applied 
to  experimental  radar  and  yawsonde  data.  The  effects  of  wind  and 
velocity  dependent  coefficients  can  readily  be  included  in  the  method. 
Position- time  radar  data  are  readily  reduced  to  aerodynamic  coefficients 
with  satisfactory  results.  Yawsonde  data  from  a  flight  test  of  the 
M549  projectile  easily  yielded  the  static  moment  coefficient.  Magnus 
and  damping  moment  coefficients  are  difficult  to  extract  from  small 
amplitude,  nutation-dominated  motion.  The  capability  of  extracting 
three-dimensional  angular  motion  from  two-dimensional  yawsonde  data  has 
been  demonstrated. 
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Schematic  of  the  155mm,  M549  Rocket-Assisted  Project 
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Small  Data  Interval  Showing  Goodness  of  Fit  with  Only  Linear  Coefficients 


Figure  5.  A  Comparison  between  Cm  Obtained  from  Yawsonde  Data  and 
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